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ABSTRACT 

It is commonly thought that stars are responsible for reionizing the Universe. However, deep searches 
for star-forming galaxies during the epoch of reionization have not yet found sufficient galaxies to 
provide the necessary ionizing flux. Fast accretion shocks associated with gravitational infall of baryons 
during the formation of galaxies have recently been proposed as an alternative method of generating 
the required ionising photons. In this Letter we assess the contribution to hydrogen reionization from 
shocked gas associated with gravitational infall. We find that shocks can ionize no more than a few 
percent of the cosmic hydrogen by z ~ 6. However, the small fraction of ionising radiation produced 
by fast accretion shocks would be significantly more biased than that associated with stars, leading to 
a modification of the luminosity weighted source clustering by ~ 10%. This modification of the bias 
may be measurable with future precision experiments utilising the redshifted 21cm line to study the 
distribution of hydrogen during the reionization era. 
Subject headings: cosmology: theory, reionization, diffuse radiation 



1. INTRODUCTION 

Star-bursting galaxies and quasars have been the lead- 
ing candidates for the sources of the UV radiation 
required to reionize th e hydrogen gas in the inter- 
galactic medium (IGM) (|Loebll2010h . The quasar pop- 
ulation is obser ved to decline quickly at z > 2.5 (e.g., 
iFan et aLll2002j ) and so it is believed that galaxies con- 
tributed the_bji!kof IJV photons that drove reion- 
ization (IMadau et al.lH99l ISrbTnovskv fc Wvithll2007t 
IBolton fc Haehneltl 12007^ 7" The observed number counts 
of high-redshift galaxy candidates (jBouwens et al.ll20TTI : 
lYan et aT1l2010[ ) discovered with the Hubble Space Tele- 
scope Ultra Deep Field (HUDF) have been used to build 
up a statistical description of star-forming activity at red- 
shift z > 7. At z ~ 7, 8.6 and 10.6 the flux limits corre- 
spond to absolute magnitudes Mn m — —18.0, —18.3 and 
— 18.6 mag, respectively. While impressively faint, these 
observations do not reach the levels corresponding to the 
faint est galaxies thought to exist at these early epochs 
(e.g.. iBarkana fc Loebl 120001 ) . and the observed stars in 
t he HUDF are i nsuffi cient to reionize the Universe. 

iTrenti et alj (|2010f ) have constructed an empirical 
model for the luminosity function and find that the ob- 
served population could have reionized the Universe if 
it extends to luminosities fainter than observed. In ad- 
dition to g alaxies, the discovery of high redshift gam ma 
ray bursts (jSalvaterra et al.l l2009: Ta nvir et al.ll2009f) can 
be used to probe t he star formati on rate up to z ~ 8 
pCistler et alJ2009(l , based on which lWvithe et all (|2010j ) 
found sufficient star-formation to achieve reionization by 
z ~ 6. 

Thus, there is evidence that stars are capable of reion- 
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izing the Universe. However, until the sources are 
identified directly, it is prudent to study alternatives. 
One such alternative is provided by cooling radiation 
associated with shock-heated gas, which can process 
gravitational energy associated with structure formation 
into an ionising radiatio n background. For example, 
iFurlanetto fc Loebl (|2004[ ) have considered the effect of 
large scale structure shocks on the IGM during reioniza- 
tion, finding that shocks heat the early IGM and produce 
a radiation background that effects molecular hydrogen 
formation and the spin temp erature of neu t ral hy drogen. 
In a complementary study, Miniati et al.1 (|2004f) calcu- 
lated the UV background from cooling radiation associ- 
ated with virialised gas that is shock heated during halo 
formation. They evaluated the ionising rate in the opti- 
cally thin IGM at z < 6, and determined that the ionising 
background at the hydrogen ionisation edge produced fol- 
lowing virialisation shocks was likely to be much smaller 
than for stars, unless supernova feedback could efficiently 
reheat the galactic gas (in which case the galaxies would 
likely produce a sig nificant contribution to the ionising 
flux from star light) . iMiniati et al.l (|2004l ) concluded that 
the hard spectrum produced may doubly reionize helium 
by z ~ 6, and also described the resulting effects on the 
thermal his tory of the IGM. 

Recently, iDopita et al.l (|2011| ) presented high resolu- 
tion simulations of the shocks associated with the super- 
sonic gravitational infall of gas onto a nascent galactic 
disk that arises when the cooling rate is much shorter 
than the free-fall time. They argued that this provided 
a source of photons that could augment the ionising flux 
from galaxies, and indeed might dominate the reioniza- 
tion of the Universe. If true this finding would have sig- 
nificant implications for studies of reionization. In par- 
ticular, the direct link between star-formation and the 
ionisation structure of the IGM would be removed, so 
that future redshifted 21cm experiments would not pro- 
vide a fruitful route to study the first stars. 

In this Letter we assess the contribution to hydro- 
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gen reionization from ionising photons produced in 
fast rad iative accretion sho cks. We note that virial- 
isation (iMiniati et al.l [2004) and fast accretion shocks 
(jDopita et al.ll2011| ) should lead to similar ionizing lu- 
minosities since the cooling radiation following the viri- 
alization shock contains energy comparable to the grav- 
itational potential energy available to drive the fast ac- 
cretion shock. In our numerical example s, we adopt the 
stand ard set of cosmological parameters (|Komatsu et al.l 
1201 11 1, with values of fi b = 0.04, Q m = 0.24 and fl A = 
0.76 for the matter, baryon, and dark energy fractional 
density respectively, h — 0.73, for the dimensionless Hub- 
ble constant, and erg = 0.82. 

2. IONIZING PHOTONS 

We begin by computing the number of ionizations per 
baryon processed through shocks. We also discuss the 
ionisation history due to stars for comparison. 

2.1. Ionizing photons based on collapsed fraction 

Dop ita et al.l ()201lD presented a fitting formula for 
the number of ionising photons that enter the IGM per 
baryon processed through a shock, which in our termi- 
nology is 



7V 7 =0.2 



iV 7 = 10 




for 100 < v < 280km/s 



for 280 < v < 400km/s (1) 



for v > 400 km/s, 



400km/s, 

where v is the velocity of the shocked gas. iDopita et al.1 
(2011) argue that v = \/2?M r where u v ir is the virial 
velocity of the halo, and we utilise this throughout the 
current work. The virial velocity is estimated from the 
halo mass using 



r (M halo ,z) -23.4 
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where £(z) is close to unity and defined as ( s 

[(ft m /Q£J(A c /i87r 2 )], n* m = [i + (n A /n m )(i + z)- 3 ]- 1 , 

A c = 18tt 2 + 82cZ - 39cP, and d = ?L z rn - 1 (see equa- 
tions 22-25 in Barkana & Loeb 2001 for more details). 

Under the assumption that gas shocks once when a 
halo forms, the number of photons produced per hydro- 
gen atom in the Univer se can be estimated from the 
rrcss fc Schechterl (fl97l mass function d n/dM (with 
modifications due to ISheth fc Tor"nTen1 IT999 N ) 
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where Mb = (fib/^m)M and pb = (^b/^m)Pm, are the 
baryonic mass inside a dark matter halo and baryonic 
mass density in the IGM respectively. Results are shown 
as a function of redshift in the upper left panel of Fig- 
ure [TJ assuming minimum halo masses corresponding to 
v v ix = 10 km s" 1 (the cooling threshold for hydrogen), 
and u v ir — 30 km s^ 1 (the Jeans threshold in an ionized 



IGM, iDiikstra et "all I2004D . The results are indepen- 
dent of this choice owing to the dominance of massive 
halos in producing radiation from fast accretion shocks 
(jDopita et al.1 l2l)llT ) . We find that only a few percent 
of the IGM is reionized by z ~ 6, and < 1% at z ~ 8, 
indicating that there are 1-2 orders of magnitude too 
few ionising photons produced in fast accretion shocks 
to reionize the Universe. 

We note that virialisation shocks do not exist in low 
mass halos due to the presence of cold flows. Given 
the halo mass dependent fraction of cold flow accre- 
tion / roiH where no shock i s pro duced (jKeres et al.l 
[20091: iFaucher-Giguere et al.l [2010), one should ex- 
clude cold mode accretion material from the calcula- 
tion of ionising luminosity associated with virialisation 
shocks (Minia ti et al.ll2004[). T his would reduce the pre- 
dictions ^irilMiniatOFalJjlOOj) by a factor of (1 — /cold), 
where / co id is weighted over the m ass and number of 
halos. However IDopita et al.l (|2011h argue that with suf- 
ficient resolution the cold flow material is found to shock 
at the intersection with the nascent galactic disk. As 
a result no correction for cold flow accretion should be 
applied to estimates of fast accretion shock produced ion- 
ising photons in equation ([3]) or later in this Letter. 

2.2. ionising photons from stars 

We compare the above result to the number of ion- 
izations obtained for stars (lower right panel of Fig- 
ure [J) . Here we utilise equation © but assume N y — 
4000(/^/ csc ) = 15, appropriate for a Salpctcr IMF with 
a fiducial value of /*/ CS c ~ 0.004, where /* and / esc 
are the star-formation efficiency and escape fraction of 
ionising photons respectively. We again show minimum 
halo masses corresponding to v v - !r = 10 km s^ 1 and 
Vv'n = 30 km s _1 , leading to significant variation in the 
total number of ionising photons produced. In difference 
to shock produced photons, these fiducial stellar popula- 
tions are easily able to reionize the Universe by z ~ 6 (as 
is well known). This result is plausible since the nuclear 
efficiency of stars is larger by many orders of magnitude 
than the efficiency of converting rest mass to radiation 
by a shock, ~ (v/c) 2 < 10~ 6 , for the shock speeds of 
interest (v < 300 km s _1 ). 

2.3. Ionizing photons based on merger rates 



The calculation in §[27T] utilises each baryon only once, 
whereas, unlike the case for stars, a baryon may be pro- 
cessed through shocks several times during the hierarchi- 
cal formation of a galaxy. We therefore recalculate the 
contribution to reionization from fast accretion shocks 
based on the merger rate of halos. Specifically, when a 
halo of mass M2 < M\ merges with a halo of mass Mi , 
we assume that the baryons contained within halo 1 are 
shocked at v = \[2v V i r {M\ + M 2 ). We assume that the 
star-formation efficiency is negligible so that all baryons 
are available to shock during each merger. The resulting 
expression for the rate of production of all baryons in 
galaxies is 
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Fig. 1. — Plots of the number of photons produced per hydrogen atom in the Universe. Upper Left: The number produced when all gas 
that enters a halo shocks once when the halo forms. Upper Right: The number produced when gas re-shocks following each merger (based 
on the merger rate of halos). Lower Left: The number produced assuming all gravitational energy of collapsed gas was released as photons 
with 13.6eV. Lower Right: The number number of ionizations obtained for stars, assuming JV 7 = 4000(/*/ eS c) = 15. Results are shown as 
a function of redshift, assuming minimum halo masses corresponding to ti v ; r = 10 km s _1 , and v v ; r = 30 km s — 1 . The grey strips represent 
uncertainty in trg. 



yielding 



dz'- 



dz 



"H 
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Here 



dzdM-2 



Mi 



is the number of mergers of a halo with 

mass between M% and M2 + dM-i that merge with a halo 
of mass M\ in a redshift interval dz. 

The assumption that baryons shock during each 
merger is optimistic. While mergers will certainly in- 
crease the virial energy per baryon, it is not clear how the 
baryons would be recycled through virialization shocks 
on multiple occasions if they cool efficiently to produce 
the ionizing radiation. If the baryons cool from the halo 
into a disk, they will not be heated to the virial tempera- 
ture of the halo (unless they have been expelled through 
some feedback process) and so should not participate in 
the halo shock of the next merger. Thus, our results 
from equation Q represent the maximum ionisation flux 
possible from the fast accretion shock mechanism sum- 
marised in equation ([T]). 

Results are shown as a function of redshift in the upper 
right panel of Figure [TJ assuming minimum halo masses 
corresponding to u v ir = 10 km s _1 and u v ir — 30 km s" 1 
(again with no discernible difference owing to the domi- 
nance of massive halos). We find that the re-processing 
of baryons through shocks in multiple mergers increases 
the ionising photon output by a factor of a few relative 



to the collapsed fraction calculation. However, we still 
find that only ~ 5 percent of the IGM is reionized by 
z 6, and ~ 1% at z ~ 8, indicating that shocks pro- 
vide insufficient ionising luminosity to reionize the IGM. 

The dif ference between our findings and the results pre- 
sented in IDopita et al.l (|2011l ) originates partly from the 
adoption of a large value of erg = 0.9 in that work, and 
partly from an error in the calculatio n method (Lawrence 
Krauss, private communication). IDopita et al.l (|201 lh 
calculate the accretion rate and corresponding ionising 
luminosity as a function of halo mass (see their Figure 
2). They then integrate over the mass- function and red- 
shift. This procedure effectively sets the rate at which 
gas in the halo doubles to be equal to the inverse of dy- 
namical time at the virial radius (which is shorter by 
an order of magnitude than the Hubble time), and so 
does not account for the duty-cycle of the shocks (which 
should be only ~ 0.1). As a result, the calculation in 
IDopita et all ((20 111 ) accretes an order of magnitude more 
gas than available per halo. 

2.4. A maximum of ionizing photons based on binding 

energy 

In the previous subsections we have estimated the num- 
ber of ionising photons available per hydrogen in the 
IGM. We next estimate the maximum number of pho- 
tons that equal the binding energy of all baryons in the 
halo. This should provide an upper limit to the number 
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of ionising photons pr oduced by shocks ( in the absence 
of significant feedback, Miniati et al. 2004), and hence an 
upper limit on the contribution of shocks to reionization. 
For this calculation, we again appeal to equation ([3]), 
setting iV 7 to be 

^ = T3^V ■ (6) 

Results are shown in the lower- left panel of Figure [T] as- 
suming minimum halo masses corresponding to v v - 1T = 
10 km s _1 and w v ir = 30 km s" 1 . We find that the total 
gravitational energy available for ionisation of hydrogen 
corresponds to less than 1 ionising photon per 3 hydro- 
gens by z ~ 6 and less than 1 ionising photon per 10 
hydrogens by z ~ 8. 

3. REIONIZATION HISTORIES 

Next we use the estimate of flux based on our merger 
calculation of ionising radiation from shocks as the 
source term in a calculation of t he reionization his- 
tory. Mirald a-Escude et al.l (]2000f ) presented a model 
which allows the calculation of an effective recombi- 
nation rate in an inhomogeneous universe by assum- 
ing a maximum overdensity (A c ) penetrated by ionizing 
photons within HII regions. The model assumes that 
reionization progresses rapidly through islands of lower 
density prior to the overlap of individual cosmological 
ionized regions. Following the overlap epoch, the re- 
maining regions of hi gh density are gradually ionized. 
iWvithe &: Loeb! (|2003l ) employed this prescription within 
a semi-analytic model of reionization, and we refer the 
reader to that paper for details of the model. Within this 
formalism, the epoch of overlap is precisely defined as the 
time when the volume fraction Q of the universe ionized 
up to an overdensity A c , reaches unity. After the overlap 
epoch, ionizing photons will experience attenuation due 
to residual overdense pockets of HI gas. The model also 
follows the mass averaged ionized fraction (Q m ). 

Figure [5] shows the resulting model for the reionization 
of the IGM and the subsequent post-overlap evolution 
due to ionising sources from fast accretion shocks (solid 
lines). Here we have used an ionising photon rate based 
on equation The case shown corresponds to a value 
for the critical overdensity prior to the overlap epoch of 
A c = 5, and both the volume averaged (dark lines) and 
mass-averaged (grey lines) ionisation fractions are shown. 
We find that shocks can reionize less than 10% of the 
IGM (by volume or mass) prior to z ~ 6, and cannot 
complete reionization until z ~ 3. For comparison we 
compute the reionization history for stars (dashed lines), 
where the ionising photon production rate is based on 
d(n~f/riii)/dz, with n 7 /nn based on equation ([3]) with 
N 7 = 15. In this model, stars complete reionization by 
z ~ 8, at which time the relative contributions from stars 
and fast accretion shocks differ by a factor in excess of 
100. 

4. IMPLICATIONS FOR THE BIAS OF IONISING SOURCES 
AND 21CM STUDIES 

In this Letter we have demonstrated that recent es- 
timates of the ionising luminosity from fast accretion 
shocks associated with galaxy formation are not suffi- 
cient to drive reionization. However the ionising pho- 
tons produced by shocks are dominated by massive ha- 



Reionization history 




Fig. 2. — Plots of the reionization history of the IGM and the sub- 
sequent post-overlap evolution due to ionising sources from merger 
induced shocks (solid lines), and stars with 7V 7 = 15 (dashed lines). 
The case shown corresponds to a value for the critical overdensity 
prior to the overlap epoch of A c = 5, and both the volume averaged 
(thick lines) and mass-averaged (thin lines) ionisation fractions are 
shown. 



los (jMiniati et alJl2004t iDopita et all 1201 lh . This is in 
contrast to the ionising radiation from stars, which is 
both predicted and observed to be dominated by low 
mass galaxies. As a result, the ionising radiation pro- 
duced in shocks is significantly more biased relative to 
the underlying large scale density of the IGM than are 
ionising photons produced in galaxies. It is easy to see 
the physics of the dominance of massive halos by not- 
ing that the collapse energy available in equation (jB]) is 
proportional to v% iT (or M^ alo ), whereas the stellar mass 
(assuming a constant mass-to-light ratio) is proportional 
to (or Mhalo). 

The ionisation structure of the IGM, particularly the 
scale of HII regions produce d is a sensitive functio n of 
the bias of ionising sources (|McQuinn et al J 120071 ). It 
is this relation between the bias of ionising sources and 
the resulting ionisation structure during reionization that 
motivates redshifted 21cm experiments with the ulti- 
mate aim of connecting galaxy p roperties to the power- 
spectrum of 21cm fluctuations (Barkanal 120091 ). Here 
we quantify the effect of fast accretion shocks on the 
bias of ionising sources. The halo bias b for a halo 
mass M at redshi f t z m ay be approximated using the 
iPress fc Schechterl (11974)) formalism mod ified to include 
non-s pherical collapse ((ShetJiet^ |200l|) . The power- 
spectrum of the space distribution of sources is propor- 
tional to b squared. The luminosity weighted bias of 
ionising radiation produced by shocks arising in mergers 
can be evaluated using the expression 
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where the bias b is evaluated at a mass Mi + Mi- The 
resulting bias is plotted as a function of redshift in the 
upper panel of Figure [3J Prior to reionization, fast ac- 
cretion shock powered ionising sources have a luminosity 
weighted bias of (6 s hock) ~ 10. 
For comparison, we calculate the luminosity weighted 
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Fig. 3. — Plots of the luminosity weighted galaxy bias. Upper 
Panel: the bias is plotted as a function of redshift due to ionising 
sources from merger induced shocks (solid lines), and stars with 
Nj = 15 (dashed lines). Lower Panel: the fractional change in 
observed clustering of ionising radiation when compared with the 
stars-only history. 



bias for stellar sources ((6 s tar)) based on the derivative of 
equation ([3]), and in analogy with equation (j4]). The re- 
sult is also plotted in the upper panel of Figure [3] Prior 
to reionization, stellar ionising sources have a smaller 
luminosity weighted bias of (6 s tar) ~ 4. Finally, we eval- 
uate the luminosity weighted bias (b) obtained when fast 
accretion shock powered ionisation sources are added to 
the stellar sources needed for reionization. We then cal- 
culate the fractional change in observed clustering of ion- 
ising radiation relative to the stars-only reionization his- 
tory [((b) 2 - (& s tar) 2 )/(&star) 2 ]- This fractional change is 
plotted in the lower panel of Figure EH which shows that 
the power-spectrum of ionising sources is increased by 
~ 10% owing to ionising radiation produced in shocks. 
This change in the clustering of ionising sources will lead 
to comparable changes in the amplitude of redshifted 
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21cm fluctuations (jWvithe fc Morales! [2007h . that will 
be detectable by pla nned low frequency radio telescopes 
(jWvithe et alJl2009h . 

5. SUMMARY 

Based on rece nt high resolution simulations 
(|Dopita et al.l 120111 ). we have quantified the contri- 
bution that gravitationally powered fast accretion 
shocks during galaxy formation can make to the reion- 
zation of hydrogen. We find that ionising radiation from 
fast accretion shocks represents a negligible contribution 
to the overall reionization history of hydrogen, leaving 
the dominant contribution to be provided by stars. 
This conclusion is consistent with expectations based 
on observations of cosmic background radiation. The 
energy released by star formation at high redshift is 
stored in the cosmic infrared b ackground with vL v ~ 10 
nW/m 2 /Str (jHauser fc Dwekll2001l ). On the other hand, 
any energy surplus from gravitational shocks is stored in 
the cosmic soft x-ray background at ~ 0.01 nW/m 2 /Str. 
This three orders of magnitude difference is suggestive 
of the relative efficiency of shocks and star formation in 
illuminating the Universe at high redshift. 

Despite their small contribution to hydrogen reioniza- 
tion, shocks may have observable conseque nces for stud- 
ies of the reionization era. As discussed by lMiniati et al.l 
(2004), the harder spectrum associated with shocks will 
lead to a modification of the thermal history. In par- 
ticular, the reionization of hydrogen by shocks would be 
accompanied by reionization of singly ionized Helium, 
thus he ating the IGM to l evels above those observed at 
z ~ 5 (jBecker et al.l 120111 ). While this likely rules out 
reionization by shocks, independently from the hydrogen 
ionisation photon budget, heating by shocks may still 
have import ant consequences for star-formation at high 
redshift (e.g. IDiikstra et alj|2004 ). 

In addition, we find that because the small contribu- 
tion from fast accretion shocks is produced in highly bi- 
ased galaxies, their presence modifies the mean clustering 
bias of the combined ionising radiation. This modifica- 
tion will likely lead to observable changes in the red- 
shifted 21cm fluctuations from neutral hydrogen during 
reionization, and so will need to be considered in anal- 
yses which aim to use precision measurements of 21cm 
fluctuations to study the properties of very high redshift 
galaxies. 
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